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INTRODUCTION 


Electromagnetic  radiation  in  the  lower  frequency  bands  emitted  in  the 
region  between  the  finitely  conducting  earth  and  the  weakly  ionized  plasma  of 
the  lower  ionosphere  is  guided  by  these  boundaries  and  propagated  to  great 
distances  around  the  earth.  The  waveguide  model  analyzes  this  propagation 
using  a  mode  series  (ref  1-4).  The  propagation  parameters  of  interest  (phase 
velocity,  attenuation  rate,  excitation  factor)  can  be  evaluated  for  each 
mode.  The  total  field  at  any  point  along  the  guide  is  then  the  vector  sum  of 
the  contributions  from  each  mode.  The  waveguide  model  has  been  employed  in  a 
series  of  computer  programs  (ref  5-7)  useful  from  elf  to  If. 

In  many  instances,  the  earth-ionosphere  waveguide  can  be  considered  to 
have  constant  propagation  properties  along  the  transmission  path.  The  mode 
sum  calculations  made  for  these  cases  are  referred  to  as  horizontally  homogen¬ 
eous.  However,  for  propagation  to  great  distances  it  is  unrealistic  to  assume 
the  waveguide  parameters  will  remain  constant  along  the  whole  length  of  the 
path.  For  example,  the  direction  and  strength  of  the  earth's  magnetic  field 
will  vary  along  the  path,  causing  changes  in  the  mode  parameters  due  to  the 
anistropy  of  the  ionosphere.  Discontinuities  can  occur  in  the  lower  wall  of 
the  waveguide  due  to  the  presence  of  ground  conductivity  changes.  The  ionos¬ 
pheric  conductivity  varies  according  to  the  time  of  day,  season,  and  the 
presence  of  the  sunrise  or  sunset  line  along  the  propagation  path.  Anomalous 
ionization  from  a  nuclear  detonation  may  also  produce  important  changes  in 
mode  parameters  along  the  propagation  path.  If  changes  in  the  waveguide 
parameters  are  small  within  a  few  radio  wavelengths,  then  the  modes  can  be 
tracked  through  the  region  of  change. 

The  computer  program  described  in  this  report  was  developed  for  calculat¬ 
ing  vlf/lf  field  strengths  in  the  earth-ionosphere  waveguide.  Allowance  is 
made  for  horizontal  inhomogeneity  in  the  direction  of  propagation.  The 
program  is  based  upon  a  slab  model  and  assumes  the  waveguide  is  homogeneous 
transverse  to  the  great  circle  path  between  a  transmitter  and  a  receiver.  For 
purposes  of  identification  it  is  called  FLD4  or  WKBFLDS. 

In  the  FLD4/WKBFLDS  program,  field  strength  calculations  can  be  generated 
for  all  electric  field  components,  at  any  receiver  height  within  the  wave¬ 
guide,  generated  by  electric  dipole  exciters  of  arbitrary  orientation  located 


/ 
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at  any  height  within  the  guide.  Transmitter  and  receiver  heights  must  be 
within  the  free-space  portion  of  the  waveguide. 

In  describing  the  WKBFLDS  computer  program,  familiarity  is  assumed  with 
references  5  and  7  which  describe  FORTRAN  programs  for  obtaining  waveguide 
mode  constants  within  the  earth- ionosphere  waveguide.  Relevant  outputs  from 
these  programs  are  the  mode  eigenangles  and  four  independent  quantities  from 
which  the  excitation  factors  may  be  determined.  Principal  outputs  of  the 
WKBFLDS  program  are  listings  of  the  mode  sum  parameters  and  the  corresponding 
plots.  The  field  strengths  are  given  in  dB  above  1  uv/m  and  the  phase  is 
given  relative  to  free  space  propagation. 

The  Summary  of  the  Equations  section  summarizes  the  relevant  formulae.  A 
description  of  the  program  input,  output,  and  operating  procedures  is  given  in 
the  Program  Execution  section.  The  appendix  contains  a  FORTRAN  listing  of  the 
computer  program. 


SUMMARY  OF  THE  EQUATIONS 


BACKGROUND 

In  the  propagation  of  vlf  and  If  terrestrial  radio  waves  to  great  dis¬ 
tances,  the  waves  are  confined  within  the  space  between  the  earth  and  the 
ionosphere.  This  space  acts  as  a  waveguide. 

The  waveguide  mode  method  described  in  references  6  and  7  obtains  the 
full-wave  solution  for  a  waveguide  that  is  characterized  by  arbitrary  electron 
and  ion  density  distributions  and  collision  frequencies  with  height  and  a 
lower  boundary  which  is  a  smooth  homogeneous  earth  with  an  adjustable  surface 
conductivity  and  dielectric  constant.  The  model  allows  for  earth  curvature, 
ionospheric  inhomogeneity,  and  anisotropy. 

The  propagation  geometry  is  shown  in  figure  1.  The  vertical  coordinate 
is  z,  the  direction  of  propagation  is  x,  and  y  is  normal  to  the  plane  of  prop¬ 
agation.  Thus,  the  fields  exhibit  no  y  dependence  and  have  a  dependence  on  x 
of  the  form  exp(-ikSmx) .  All  field  quantities  are  assumed  to  have  an  exp(iwt) 
dependence  where  ui  is  the  angular  frequency.  The  dipole  source  for  the  fields 
is  denoted  in  figure  1  by  The  dipole  is  oriented  at  an  inclination  angle  Y 
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measured  from  the  vertical,  and  azimuthal  angle  0  measured  from  the  x-axis. 
For  a  horizontal  dipole  Y  =  90° ,  and  0=0°  represents  end-fire  launching. 

The  energy  within  the  waveguide  is  considered  to  be  partitioned  among  a 
series  of  modes  called  the  eigenangles  (or  "modes")  and  designated  by  the 
complex  angle,  0,  obtained  by  solving  the  determinantal  equation: 

F( 0)  =  | R( 0)  1(0)  -  1 |  =  0  (1) 


where 


R  (0)  = 

and 

R(0)  = 


1|V0) 

iV0) 

R  (0) 

II  1 

11  / 

R  (0) 

°  \ 

II  U 

1R1(G)  / 


(2) 


(3) 


are  the  conqplex  ionospheric  reflection  coefficient  matrices  looking  up  into 
the  ionosphere  and  down  towards  the  ground,  respectively,  from  a  height  "d". 

Once  equation  1  is  solved  for  as  many  modes,  ©m,  as  needed,  the  mode 
parameters,  phase  velocity,  attenuation  rate,  and  excitation  factors  can  be 
computed.  Let  Sm  be  the  sine  of  the  eigenangle  0m  and  let  the  individual 
components  be  Sr  and  so  that  Sm  =  Sr  +  i  The  phase  velocity  is  given 

by 


c 


where  c  is  the  vacuum  speed  of  light. 
The  attenuation  rate  is  given  by 


a  =  -8.6859k 


(in  dB/unit  distance) 


where  k  is  the  free-space  wave  number. 


(4) 


(5) 
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The  excitation  factor  formulae,  as  given  in  reference  1,  are  summarized 
in  table  1.  The  headings  apply  to  the  electric  field  components  (E z,  Ey,  Ex> 
excited  by  a  vertical  dipole  (z),  a  horizontal  dipole  broadside  (y),  and  a 
horizontal  dipole  end-on  (x). 


Field 

Component 


Exciter 


s5/2  nyy  (i-^y 

F’(0m)  11*11  °11 


3/2  llR1(1+(|R||)(1+iRi)  s3/2  (  1+  , R  „  ) 2  C  1- jR^) 


F'  (G  ) 
m 


12 


F’  (0  ) 
m 


n*n  Dn 


s3/2  sV2 


y  f'  (0  )  d 

m  12 


F’(Q  )  .R.  D 

m  i  1  22 


s1/2  ^,(1+^(1+ |R|) 


F'  (0  ) 
m 


12 


^  .2 


s3/2  <1+„V  ^llV  s1/2  l|Ri(  1+  (|R  „ )  (  1+iRi) 


F'(G  )  ..R,,  D 

m  II  II  11 


E"  (0  ) 
m 


12 


.1/2  (  1+  l|R u >  <  1_  iRlxRi) 


F'(©  >  ,R,  D 

m  II  II  1  1 


Table  1 .  Excitation  factors 


In  table  1  the  following  terms  are  used: 


r'<v=H 


q=0 


'r'V3‘<  -r2>’/2  wl 


.  no  ,ak,1/3,_  2  2 . 1/2 


F2  =  H1«*o>  ^T1?  %  -  S  >  V^o* 


,  _  (W  . ,  .ii/V  sV"  vvi 


(6) 


(7) 


(8) 


(9) 
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(  10) 


f  -  h;<q  ,  -  i(^,'/3,n2  -  s2,’/2  v,  , 

4  1  o  2  g  1  o 


2  ~2/3  2  2 
q  =  (  — )  (1  +  -  z  -  S) 

z  a_k  a 


(11) 


2/3 

H.(q)  =  h'(q)  (Jj^)  h  (q)  ;  j  =  1,2 


(  12) 


2  ,  ^  2 
n  =  1  +  — z 
z  a 


(13) 


n  =  - - l  - 

g  e 

o  o 


(14) 


D11  ^iN^d*  +  F2h2(qd^ 

°12  =  KW  +  F2h2(qd^  tP3h1(V  +  F4h2(qd)} 

D22  =  iF3h1(qd}  +  F4h2(qd^2 


(15) 

(16) 

(17) 


where  e/e0  =  dielectric  constant  of  the  ground  and  o  =  the  ground  conductiv¬ 
ity. 


The  subscripts  on  q  and  n  represent  the  value  of  z  at  which  these  quantities 
are  evaluated.  The  functions  h^  and  h2  are  modified  Hankel  functions  of  order 
1/3  as  defined  in  reference  8.  The  primes  on  these  quantities  denote  deriva¬ 
tives  with  respect  to  the  argument. 

The  modal  excitation  factor  and  height  gain  functions  are  needed  in  com¬ 
puting  electric  field  strengths.  The  excitation  factors  of  table  1  must  be 
supplemented  with  the  height  gain  functions  which  are  given  in  reference  1  as: 


f  (z)  =  exp  (-)  {f  h  (q)  +  F  h  (q)} 
I  ail  z  z 


(18) 
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(  19) 


f2(Z)  =  ^  ^3h  1 ( ^  >  +  F4h2(q)^/^F3h1(qo)  +  F4h2<qo^ 


f3(z) 


Jl_ 

ik 


df , 


dz 


(20) 


where  is  the  height  gain  for  the  vertical  electric  field  E  ,  f2  is  the 

height  gain  for  the  horizontal  electric  field  component  E  ,  and  f,  is  the 

y  ^ 

height  gain  for  the  horizontal  electric  field  component  E^. 

The  basic  mode  sum  equation  is  given  by 

£  E 

—  (21) 

[a  sin(x/a)] 

where  E  =  total  electric  field  in  volts/metre 
Em  =  complex  amplitude  of  each  mode 
p  =  power  radiated 
H  =  free-space  impedance 


F  =  H  r  1 1/2  k 
4TT  ''10k'1  2 


If  the  waveguide  is  homogeneous  in  the  direction  of  propagation,  then  E  is 

ra 

given  by 

E  =  A  exp  {-ik  S  x}  (22) 

m  m  1  m  j 

where  A^  is  the  complex  amplitude  of  the  mth  mode  to  be  described  later.  If 
the  parameters  Am  and  Sm  are  slowly  varying  functions  of  distance,  then  the 
so-called  WKB  form  of  is  given  by  references  1  and  9  as 


where 


m  ,  v ,  1/2 


An  =  UAzr(0)  Azr(x)!  gz(zt>  cosY 


+  [Am  (0)  Am  (x)] 1/2  g  ( z  )  sinY  sin6 
yr  yr  y  t 

+  tVjjr(°)  Axr  (x)!1//2  9X(V  sinY  cose)  9r  ( zr ) 


(24) 


and  r  denotes  the  orientation  of  the  received  electric  field  component.  If 
equation  24  is  used,  then  equation  26  reduces  to 


A  -  {Am  g  (z  )  cosY 
m  1  zr  ^z  t 


U( 

+  A  g  (z  )  sinY  sin0 
yr  y  t 


+  Axr  9x(zt}  sinY  cos6l  9r(zr) 


(25) 


KB  IMPLEMENTATION 

The  procedure  implemented  in  the  WKBFLDS  conputer  program  is  to  segment 
the  earth-ionosphere  waveguide  into  M  cascaded  uniform  waveguides  as  shown  in 
figure  2.  These  segments  (or  slabs)  are  numbered  from  left  to  right.  Slab  1 
contains  the  transmitter.  The  distance  from  the  transmitter  to  the  beginning 

of  each  slab  is  denoted  by  p^<  The  waveguide  parameters  associated  with  each 
slab  are  determined  by  the  user  based  on  consideration  of  the  variations  of 
the  path  conditions  (geomagnetic  field,  ground  conductivity,  and  the  ionos¬ 
pheric  profiles  of  electron  density,  ion  density,  and  collision  frequency). 

The  mode  generating  programs  (ref  6  or  7)  are  run  to  obtain  as  many  modes 
as  desired  in  each  slab.  The  WKB  model  requires  an  equal  number  of  modes  in 
each  slab  and  the  modes  ntist  be  traced  from  one  slab  to  the  next  in  order  to 
preserve  mode  numbers.  At  each  slab  boundary  energy  is  transferred  from  a 
given  mode  into  the  corresponding  mode  in  the  next  slab.  The  values  of  A^  an(j 
at  the  beginning  of  each  segment  are  defined  to  be  those  obtained  from  the 
full-wave  programs.  The  values  of  each  quantity  within  each  segment  are 
obtained  by  linear  interpolation  between  the  values  at  the  segment  bound- 


Figure  2.  Diagram  illustrating  propagation  path  segmentati 


aries.  Care  is  taken  to  insure  continuity  of  the  phase  of  Am  and  Sm  from  the 
transmitter  to  the  receiver. 


PROGRAM  EXECUTION 


OPTIONS 


The  program  has  three  output  options.  For  compatibility  with  other  pro¬ 
grams  these  options  are  numbered  1,  2,  and  4.  Option  1  produces  fields  as 
functions  of  distance  parametric  in  transmitter  altitude  and  orientation  for  a 
fixed  receiver  altitude  and  orientation.  The  program  generates  501  points 
between  the  transmitter  (distance  =  0)  and  the  end  of  the  path  (distance  = 
DMAX ,  described  below).  The  calculated  fields  and  appropriate  constants  may 
be  written  to  logical  file  2  for  processing  by  other  programs.  The  calculated 
amplitude  and  relative  phase  may  be  printed  and/or  plotted.  The  vertical  and 
horizontal  (distance)  dimensions  and  scales  of  the  plots  may  be  arbitrarily 
specified.  Each  set  of  distance  calculations  can  be  made  for  up  to  20  trans¬ 
mitter  orientations. 

Option  2  produces  fields  as  functions  of  a  transmitter's  position  in  an 
orbit  for  a  fixed  transmitter  altitude  and  antenna  inclination  and  for  a  fixed 
receiver  altitude  and  orientation.  The  program  generates  73  points  between  9 
=  0°  and  0  =  360°.  Up  to  three  different  orbit  configurations  may  be  ob¬ 
tained:  a  dipole  rotating  in  a  counterclockwise  direction,  and  clockwise  and 
counterclockwise  orbiting  at  a  specified  radius.  Clockwise  is  determined  by 
looking  down  the  z-axis  towards  the  x-y  plane.  Each  set  of  orbit  calculations 
is  made  at  as  many  as  20  different  distances  from  the  transmitter.  The  calcu¬ 
lated  amplitude  and  relative  phase  may  be  printed.  The  results  are  always 
plotted  using  a  fixed  vertical  size  and  scale  and  a  user  specified  horizontal 
size. 

Option  4  produces  fields  as  functions  of  distance  for  fixed  values  of  the 
transmitter  and  receiver  altitudes  and  orientations.  The  data  sets  are  varied 
so  as  to  produce  comparison  plots.  The  printed  and  plotted  outputs  are  all 
specified  as  for  option  1.  The  procedure  for  storing  the  calculated  fields  on 
an  external  device  is  also  provided  as  for  option  1. 
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CONTROL  CARDS 


The  flow  of  program  execution  is  controlled  by  a  set  of  cards  which 
define  the  type  of  data  being  input  and  allow  standard  defaults  to  apply.  All 
80  columns  of  these  control  cards  are  read  and  printed.  However,  only  the 
first  four  columns  are  examined  by  the  program.  Thus  the  input  card  and  the 
printout  can  contain  additional  comments.  For  example,  "NAMELIST  DATA  FOR  AN 
ELEVATED  TRANSMITTER"  can  be  used  in  place  of  the  minimum  requirement  of 
"NAME".  The  control  cards  are  described  below  and  a  sample  set  is  shown  in 
figure  3. 


NAME 

I  OPT 
I  COMP 


TALT 

RALT 

INCL 

THETA 

NRA 

DIST 

NRD 

NRDATA 

RADIUS 

POWER 


Signals  that  NAMELIST  data  follow.  The  NAMELIST  name  is  DATUM. 
The  NAMELIST  variables  are  described  below. 

Option  number  -  DEFAULT  =  1 
Index  of  received  electric  field  component 
=  1  gives  the  vertical  field  (z)  -  DEFAULT 
=  2  gives  the  horizontal  transverse  field  (y) 

=  3  gives  the  horizontal  longitudinal  field  (x) 

Transmitter  altitude  in  km  (ie,  zt>  -  DEFAULT  =0.0 
Receiver  altitude  in  km  (ie,  z )  -  DEFAULT  =0.0 

Transmitting  antenna  inclination  in  degrees  (ie,  Y).  This  is  an 
array  of  20  elements  -  DEFAULT  =0.0 

Transmitting  antenna  azimuth  in  degrees  (ie,  9).  This  is  an 
array  of  20  elements  -  DEFAULT  =0.0 

Number  of  pairs  of  INCL  and  THETA  to  be  used.  This  applies  only 
to  option  1  -  DEFAULT  =  1 

Distance  between  an  orbiting  transmitter  and  a  receiver  in  Mm. 
This  is  an  array  of  20  elements  -  DEFAULT  =1.0 

Number  of  values  of  DIST  to  be  used.  This  applies  only  to 
option  2  -  DEFAULT  =  1 

Number  of  input  data  sets  to  process  for  option  4  -  DEFAULT  =  1 
Radius  of  the  orbit  in  Mm  used  in  option  2  -  DEFAULT  =0.0 
Total  radiated  power  in  kW  -  DEFAULT  =1.0 
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NPRINT  Printout  control  flag 


NAPLOT 


NPDIFF 

NRCURV 

SIZEX 

SIZEY 

SIZEX 1 
SIZEY 1 

SIZEX2 

SIZEY2 

AMPMAX 

AMP MIN 
AMPINC 
PHSMAX 

PHSMIN 

PHSINC 

DMAX 

DMIN 


=  0  generates  a  minimum  of  print 

=  1  generates  full  print  (described  in  example  data  section)  - 
DEFAULT 

Amplitude  plot  flag 
=  0  deletes  plot 
=  1  generates  plot  -  DEFAULT 
=  0  deletes  plot  -  DEFAULT 
=  1  generates  plot 
Phase  difference  plot  flag 
=  0  deletes  plot  -  DEFAULT 
=  1  generates  plot 

Number  of  curves  to  be  plotted  per  graph  for  option  1  and  4  - 
DEFAULT  =  4 

Length  of  the  horizontal  plot  axis  fo*-  .urrent  option  in  inches 
(for  all  options)  -  DEFAULT  =  10.0 

Length  of  the  vertical  axis  for  current  option  in  inches  (for 
options  1  and  4  only)  -  DEFAULT  =9.0 

Length  of  axis  for  option  1  and  4  only  -  DEFAULT  =  10.0,  8.0 

Lengths  of  axes  for  option  2  only  -  DEFAULT  =  1.0,  4.0  (note 
that  SIZEY2  is  included  for  completeness  only,  the  value  cannot 
be  changed) 

Maximum  value  of  amplitude  to  be  plotted  in  dB  above  1  PV/m  for 
options  1  and  4  -  DEFAULT  =  70 .0 

Minimum  value  of  amplitude  to  be  plotted  -  DEFAULT  =  -10.0 
Tic  mark  intervals  for  amplitude  axis  in  dB  -  DEFAULT  =  10.0 
Maximum  positive  phase  excursion  to  be  plotted  in  degrees  for 
options  1  and  4  -  DEFAULT  =  360 .0 

Maximum  negative  phase  excursion  to  be  plotted  -  DEFAULT  = 
-360  .0 

Tic  mark  intervals  for  phase  axis  in  degrees  -  DEFAULT  =  90.0 
Maximum  distance  of  plot  axis  in  Mm  for  options  1  and  4  - 
DEFAULT  =  10.0 

Minimum  distance  in  Mm  -  DEFAULT  =  0 .0 

Tic  mark  intervals  for  distance  axis  in  ton  -  DEFAULT  =  1.0 


XINC 


TOT APE 

TLONG 

TLAT 

RBEAR 

After  reading 

DATA 


Integer  flag  for  writing  calculated  fields  to  logical  unit  2 
=  0  deletes  output  -  DEFAULT 
=  1  generates  output 

Transmitter  longitude  (written  to  logical  unit  2  when  TOTAPE  = 
1 )  -  DEFAULT  =  0 .0 

Transmitter  latitude  (written  to  logical  unit  2  when  TOTAPE  =  1) 
-  DEFAULT  =0.0 

Geographic  bearing  from  transmitter  to  receiver  (written  to  unit 
2  when  TOTAPE  =  1 )  -  DEFAULT  =0.0 

the  NAMELIST  the  program  reads  the  next  control  card. 

Signals  that  the  propagation  path  data  follow.  The  format  of 
these  data  is  that  which  is  produced  by  reference  7.  The  first 
card  contains  the  data  set  identification.  This  card  is  used  to 
supply  a  literal  path  description  such  as  "HAWAII  TO  WAKE". 

This  is  followed  by  sets  of  mode  constants,  one  for  each  path 
segment.  The  first  card  for  each  segment  contains  the  starting 
distance  (P),  frequency  (f),  magnetic  azimuth,  codip  and 

intensity,  ground  conductivity  (0),  dielectric  constant  (e/e0), 
and  the  nominal  height  of  the  free-space  portion  of  the 
waveguide  (h).  Of  these  parameters,  P,  f,  o,  and  c/£0  must 

appear.  The  magnetic  parameters  and  h  are  included  for 

reference.  This  card  is  followed  by  the  mode  constants,  one 
pair  for  each  mode,  containing  the  following  information: 

1  e  t1  t2 

2  Q  T3  T4 

where  1  and  2  are  sequencing  indices,  9  is  the  complex  eigen- 
angle  at  the  ground,  and  Ts  are  complex  constants  described 
be  low. 


T  = 

1 


1/2-2  - 
S  (1+  R)  (1-  R  R  ) 
II  IT  1111 


F,((Y  „R,i  D,i 
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In  terms  of  the  Ts  the  elements  of 


(A  ) 

ij 
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-T  T  S  T 

3  4  2 


table  1  are  given  as: 


T  S 
1 


\ 


-T  T 
3  4 


(30) 


1  -T  S  T  -T  / 

\  1  3  1  / 

The  list  of  inodes  for  each  segment  is  terminated  by  a  blank 
card.  A  maximum  of  30  modes  will  be  used  although  there  is  no 
maximum  number  which  the  data  deck  may  contain.  The  list  of 
segments  is  terminated  by  a  card  with  p  =  40 .  Parameters  for 
each  segment  are  stored  sequentially  on  logical  file  3.  The 
number  of  path  segments  is  limited  by  the  space  allocated  to 
this  file.  Each  segment  requires  1445  words  of  storage.  After 
reading  these  data  the  program  returns  for  another  control  card 
if  the  option  is  1  or  2.  If  the  option  is  4,  then  the  program 
generates  the  required  calculations.  Upon  completion  of  these 
calculations,  the  program  returns  to  read  the  next  data  set 
starting  with  the  data  set  identification  card.  This  cycle 
continues  until  NRDATA  sets  have  been  processed. 


START  Signals  that  all  input  is  complete  and  initiates  execution  for 
options  1  and  2. 

A  sample  of  how  the  control  cards  may  be  used  is  shown  in  figure  3. 
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SPECIAL  CONSIDERATIONS 


The  maximum  value  of  NRCUFV  is  4  and  the  maximum  number  of  modes  is  30. 
If  either  of  these  values  is  exceeded,  a  message  is  printed  and  execution 
continues  using  the  maximum  value. 

Option  1:  Calculations  all  begin  at  the  transmitter  and  end  at  DMAX . 

DMAX,  DMIN,  and  XINC  are  used  for  the  horizontal  axis  scaling. 

Option  2:  Scaling  for  amplitude  and  phase  is  done  automatically.  The 

range  on  the  horizontal  axis  is  0°  to  360°  in  orbit  angle.  If  RADIUS  =  0.0, 

only  the  rotating  dipole  calculations  are  made  for  counterclockwise  rotation. 

Option  4:  NAMELIST  input  (specifying  option  4  and  NRDATA)  must  precede 

the  DATA  card.  Each  set  of  mode  constants  must  be  preceded  by  an  identifica¬ 

tion  card  and  not  by  a  DATA  card.  Calculations  all  begin  at  transmitter  and 
end  at  DMAX  using  the  first  value  of  INCL  and  THETA  in  each  list. 

TOTAPE  =  1  :  Unit  2  receives  the  500  values  of  the  mode  sum  as  single 

precision  complex  values,  freq,  TLONG ,  TLAT,  RBEAR,  POWER,  INCL,  THETA,  TALT, 
RALT,  DMIN,  and  DMAX  (in  that  order).  The  output  is  written  in  unformatted 
form. 

The  procedure  for  using  homogeneous  sections  along  a  path  is  to  include  a 
single  card  in  the  form  'R  XX .XXX'  after  the  appropriate  segment  data.  This 
indicates  that  the  previous  set  of  mode  constants  are  to  be  used  at  the  dis¬ 
tance  given  by  XX. XXX  on  the  card.  As  a  special  case,  when  running  the  pro¬ 
gram  for  the  case  of  a  horizontally  homogeneous  waveguide,  mode  constai.t  data 
are  input  at  'R  0.0'.  This  is  followed  by  a  single  card  in  the 
format' R  XX. XXX'  where  XX. XXX  is  the  value  of  DMAX.  Furthermore,  the  program 
allows  for  an  instantaneous  change  in  path  mode  constants.  This  is 
accomplished  by  using  the  same  value  of  P  for  two  consecutive  sets  of  mode 
constants.  This  is  used  for  minimizing  the  discrepancies  generated  by  the 
program  when  the  WKB  assumption  is  violated,  such  as  for  large  changes  in 
ground  conductivity.  An  example  of  the  horizontally  homogeneous  sections  is 
shown  in  the  sample  data  set  of  figure  3. 
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1 


1 

NAME 

2 

SDATUM 

3 

10PT«1 , 

4 

JNCI«0. ,90 

.,  THETA 

l.{ 

).  ,90.  ,  NRA «2 , 

S 

SIZEX«5. , 

SIZEY*4 . 

1 

6 

SEND 

T 

OATA 

8 

TEST  OATA 

9 

R  0.0  F 

10.0 

A 

0.0  C  0.0  M  0.0 

S  10.0  E  0.0 

10 

1  90.0 

0.0 

1 

0.0  -4.70745100-002 

0.0  -4.70745100-002 

11 

2  90.0 

0.0 

1 

0.0  -2.35372600-002 

1.0  0.0 

12 

1  81 .93069 

0.0 

1 

1 .38439500-002-1 .90354300-002 

1 . 38439500-002-1 . 90354300-002 

13 

2  81 .93069 

0.0 

1 

6.92197500-003-9.51771500-003 

1.0  0.0 

14 

IS 

R  5.0 

16 

R  5.0  F 

10.0 

A 

0.0  C  0.0  M  0.0 

S  10.0  E  0-0 

IT 

1  90.0 

0.0 

1 

0.0  -4.70745100-003 

0.0  -4.70745100-003 

18 

2  90.0 

0.0 

1 

0.0  -2.35372600-003 

1.0  0.0 

19 

1  81 .93069 

0.0 

1 

1 .38439500-003-1 .90354300-003 

1 . 38439500-003-1 . 90354300-003 

20 

2  81 .93069 

0.0 

1 

6.92197500-004-9.51771500-004 

1.0  0.0 

21 

22 

R  10.0 

23 

R  40. 

24 

START 

25 

NAME 

26 

ADATUM 

2T 

ICOMP  =  2 , 

28 

NPPlOT=1 , 

29 

I NCL  =  0 . ,90 

. ,90. , THETA- 0 . ,0. ,90. ,  NRA«3 , NRCURV-3 , 

30 

SEND 

31 

START 

32 

NAME 

33 

SOATUM 

34 

ICOMP* 1 , 

35 

IQPT*  2 , 

36 

INCL  =  45 . , 

37 

DIST=4.9,5 

. 1 ,  NRD«2, 

38 

SEND 

39 

START 

40 

NAME 

41 

4DATUM  IOPT«4,NRDATA»2.NPPIOT*0,4END 

42 

OATA 

43 

TEST  OATA 

44 

R  0.0  F 

10.0 

A 

0.0  C  0.0  M  0.0 

S  10.0  E  C.O 

45 

1  90.0 

0.0 

1 

0.0  -4.70745100-002 

0.0  -4.70745100-002 

46 

2  90.0 

0.0 

1 

0.0  -2.35372600-002 

1.0  0.0 

47 

1  81.93069 

0.0 

1 

1.38439500-002-1 .90354300-002 

1.38439500-002-1 .90354300-002 

48 

2  81 .93069 

0.0 

1 

6.92197500-003-9.51771500-003 

1.0  0.0 

49 

50 

R  5.0 

51 

R  5.0  F 

10.0 

A 

0.0  C  0.0  M  0.0 

S  10.0  E  o.o 

52 

1  90.0 

0.0 

1 

0.0  -4.70745100-003 

0.0  -4.70745100-003 

53 

2  90.0 

0.0 

1 

0.0  -2.35372600-003 

1.0  0.0 

54 

1  81.93069 

0.0 

1 

1.38439500-003-1 .90354300-003 

1 .38439500-003-1 .90354300-003 

55 

2  81 .93069 

0.0 

1 

6.92197500-004-9.51771500-004 

1.0  0.0 

58 

57 

R  10.0 

50 

R  40. 

59 

TEST  OATA  WITH  Alt 

AMPLITUDES  REDUCED  BY  FACTOR  UF 

10 

60 

R  0.0  F 

10.0 

A 

0.0  C  0.0  M  0.0 

S  10.0  E  o.o 

61 

1  90.0 

0.0 

1 

0.0  -4.70745100-003 

0.0  -4.70745100-003 

62 

2  90.0 

0.0 

1 

0.0  -2.35372600-003 

1.0  0-0 

63 

1  81.93069 

0.0 

1 

1 .38439500-003-1 .90354300-003 

1 .38439500-003-1  .90354300-003 

64 

2  81 .93069 

0.0 

1 

6.92197500-004-9.51771500-004 

1.0  0.0 

65 

68 

R  5.0 

67 

R  5.0  F 

10.0 

A 

0.0  C  0.0  M  o.o 

S  10.0  E  0.0 

68 

1  90.0 

0.0 

1 

0.0  -4.70745100-004 

0.0  -4.70745100-004 

69 

2  90.0 

0.0 

1 

0.0  -2.35372600-004 

1.0  0.0 

70 

1  81.93069 

0.0 

1 

1.38439500-004-1 .90354300-004 

1 .38439500-004-1 .90354300-004 

71 

2  81 .93069 

0.0 

1 

6.92197500-005-9.51771500-005 

1.0  0.0 

72 

73 

R  10.0 

74 

R  40. 

t _  Card  column  1 


Figure  3.  Sample  data  deck 
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In  order  to  facilitate  mode  matching  across  segment  boundaries,  a  missing 
mode  can  be  indicated  by  a  pair  of  filler  cards  inserted  in  the  deck  where  the 
mode  would  have  appeared.  This  pair  of  cards  must  contain  a  non-zero  eigen- 
angle  value.  Typically  the  eigenangle  is  the  value  from  some  previous  segment 
which  had  the  mode.  The  rest  of  each  of  the  cards  is  blank. 


EXAMPLE  DATA 


The  data  listed  in  figure  3  exercise  the  program  in  most  of  its  options. 
The  test  data  are  fictitious.  There  are  only  two  modes;  their  attenuation 
rate,  phase  velocity  relative  to  the  speed  of  light,  and  complex  excitation 
factor  A  are  given  in  table  2. 


Mode 

1 

2 


Attenuation  Rate 
dB/Mm 
0 
0 


Phase  velocity/c 

1 .00 
1  .01 


dB 

0 


-6 


<pzz  in  degrees 

0 

36 


Table  2.  Test  data  mode  parameters 


These  modes  produce  interference  minima  at  3  Mm  intervals.  The  values  of  the 

Ts  of  equations  28-31  are  defined  such  that  A  =A  =  A  =  A  =  A  and 
^  zz  yy  xx  xz  zx 

A  =  A  =  A  =A  =  Vo  A  for  each  mode.  Thus,  a  horizontal  broadside 
yz  zy  xy  yx  zz 

antenna  produces  fields  which  are  the  same  amplitude  as  produced  by  a  vertical 
antenna.  In  figure  3,  lines  1-6  are  the  NAMELIST  input  data  specifying  two 
transmitting  antenna  orientations:  vertical  and  horizontal  broadside.  Lines 
7-23  are  the  propagation  path  constants.  Lines  15  and  22  illustrate  the  use 
of  horizontally  homogeneous  segments.  The  first  two  pages  of  output  for  the 
calculations  specified  in  lines  1-24  are  shown  in  figures  4  and  5.  In  figure 
4  we  see  the  list  of  segment  parameter  cards  after  the  control  card  DATA. 
After  the  control  card  START  we  see  a  listing  of  the  mode  parameters  for  each 


20 


8 

-* 

O 

O 

O 

O 

0 

-O  O 

0 

0 

0 

O  — 

0  - 

0  - 

0  • 

0 

0 

0 

0 

0 

0 

<  ID 

v  in 

v  in 

v  in 

« 

0 

0 

0 

d  0 

d  0 

d  0 

d  n 

0 

0 

0 

UJ  *  * 

Ui  •  • 

Ui  *  • 

Ui  •  • 

0  0 

0  0 

0  0 

in  — 

in  — 

in  — 

in  — 

X  d 

0  0 

0  0 

•  0 

f-  GO 

r-  cd 

OD  0> 

m  0) 

UiO  * 

•  0 

•  0 

0 

1  1 

1  1 

1  1 

1  1 

NON 

0 

0 

-  0 

—  40 

-0 

-0 

0 

WOO  1 

0 

0 

0 

■  *  0  4 

0 

0 

0 

—  O  O  #- 

0 

0 

0 

Ui  O  O  O  O 

• 

• 

0  • 

—  d 

—  d 

—  d 

—  d 

Q-  O  O  O  -# 

0 

0 

0 

t-  OD 

f*  03 

r*»  cs 

r*  0 

<  4  0  0  a 

0  • 

0  - 

0  - 

m  •  - 

CD  *  * 

0  .  . 

0 

•-  0  0  0  a. 

0 

0 

0 

a>  <x 

CD  < 

GD  < 

CD  < 

0  0  —  0  z 

0 

0 

0 

1  *- 

1  — 

—  d 

—  d 

#-o*o* 

0 

0 

0 

1 

1 

1  1 

(  1 

-  0  1  —  — 

0 

0 

0 

—  0  *  - 

0  0 

0  0 

0  0 

0  0 

0  0 

0  0 

•  0 

0  0  z  • 

•  0 

•  0 

0 

♦  -• 

O  •  2  * 

o  a  < 
o  ■  55 
o  <  o 
o  >  •  - 

O  UJ  d  d 


O  O 

•  o  . 

o 

-  o  * 

o 
o 
o 
o 

o  o  o  * 

o  •  o 

o  O  < 

O  -  O  #- 
o  O  < 

o  o  o 

o  o  cr 

o  o  z 

.  o  •  - 

O  d 

.  o  * 

o 
o 
o 
o 


0  x  <  a.  • 

00  00 

0  0 

Ui 

OUJ  •  • 

0  0  0  0 

•  0 

0 

oiM  —  r>#— 

•  0  •  0 

0 

O 

0  ~  0  0  z 

0  O 

►  0 

2 

•  in  O  O  — 

*  0  -0 

0 

-  4  4  a 

0  *  0 

O  N 

0 

*  0  0  a 

0  0 

O 

• 

0  0  z 

00  0 

O  < 

#-  0  0  - 

-  cr 

0  -  cc 

o  o  <r 
<000 
cr  o  o  o 
•  o  o  U> 
o  <  m 
o  •  • 

O  I 
OHM 
O  O 

o  •  n  z 

o  >  •- 

O  »  ui  2  ■ 

O  NU) 


o  z 
•  -o  -  - 
o 
o 
o 
o 

o  o  o  o  o 
00  •  o  o 
*00  00 

00  -00 
-00  00 

00  00 

00  00 


d  0 

d  0 

y 

1  | 

1  1 

u 

0 

T5 

<D 

CD 

0  r- 

t — < 

a 

in  d 

0  d 

n  r- 

L.r 

0  r-* 

in 

d  — 

a 

d  — 

1  l 

—4 

1  1 

u. 

W-i 

0  in 

0  in 

0 

0 

O  0  < 

d  0 

Z 

d  0 

ry 

in  — 

ui 

0  — 

0  0 

0  0 

r—i 

1  — 

1  — 

1 

• 

l/> 

O 

d  n 

d  0 

c 

- 

v  — 

•H 

-  in 

UJ 

—  in 

r-i 

O 

•  • 

n  d 

1 

in 

0 

0  d 

>1 

£> 

d  0 

< 

d  0 

T3 

0 

r-  0 

cr 

r-  0 

03 

•  • 

0 

•  • 

U 

CD  < 

0  <• 

O 

d  0 

d  0 

T5 

0  cc  -•  z  > 

0  0 

0  0 

0  <  in  a  cr 

•  • 

•  -  0 

0  Ui  -  -3 

0  0 

0 

•  0  —  0  0 

0  -  -  0 

-  -0 

a  0  0  cr 

0  0 

0 

■  -0  0  z 

-0  0 

0 

4  4- 

0  0 

0 

»-  0  0  — 

0  0 

0 

-j  000  ©O' 

<  000  o  o  < 

#-000  +  0  *  o  < 

•OOOOO  O' 
*-00000  -  o  < 

OOOOO  O  I 
00000  O' 
O  •  -00  O' 

o  00 

o  n  ■  o  O  O 

*■  «o  * 

x  •  o 

»  ►<  •  o 

Ui  2  II  o 

►-  rsj  in  is)  O 

«  <*-*10300' 
-10  Q.  Z  o  O  ' 

a  ♦-  -  -  O  *0' 

Z  -*-  «  X  <  O' 

o  00  -  a  -  o 

(J  00  -—O' 

M  4  4  OOOl 
-  00  O  o • 

—  000  ♦  • 

0000  o 

0000  o  - 

0000  o 

0000  o 

0000  o 

O  —  —  O  o 

0**0  O  O  ' 

•  O  —O' 


>000000 
>000  •  o  o 

>  *00  00 

>  00  -00 

>  -  o  o  d  o  o 

>  d  o  o  o  o  o 
>000000 
•  000  +  00 
*  4-  •  •  O  •  • 

o  o 

* o  -  •©  -  - 
o  o 

o  o 

o  o 

o  o 

0  0  0  0)00 
0)0  0  •  o  o 

•  00  00 

00  -00 

-00  00 

00  00 

00  00 

00  00 

o  o 

>  o  -  -  o  -  • 

o  o 

o  o 

o  o 

o  o 

000000 

000  *00  < 

•  00  00  ♦- 

00*00  < 

■00  00  o 

o  o  <  o  o 

—I  O  O  f-  o  o  o  <  *- 

OOOUiOOZ  ♦-  M 

2  *  ■  I  •  *  ui  <  Ui 

-I  #-  4*  O  *- 


in 

in 

in 

in 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

• 

• 

• 

* 

2 

2 

2 

2 

O 

0 

0 

0 

0 

O 

0 

z 

O 

0 

O 

0 

cc 

O 

0 

in 

3 

O 

«•* 

O 

u 

> 

(J 

0 

UJ 

O 

0 

CC 

O 

0 

O 

0 

a. 

0 

0 

O 

0 

0 

0 

• 

• 

2 

• 

» 

0 

cr 

u. 

< 

< 

< 

< 

< 

#- 

O 

0 

< 

0 

0 

O 

0 

O 

0 

0 

O 

0 

0 

0 

O 

0 

O 

0 

0 

« 

• 

z 

• 

• 

O 

0 

m 

«- 

3 

u. 

u. 

u. 

u. 

O 

0 

0 

0 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

*" 

K 

c 

K 

K 

0  0 

O  0 

O  0 

0  d 

0  d 

0  d 

O  0 

O  0 

O  0 

3 

_J  *  * 

_J  .  . 

_J  -  • 

a 

< 

< 

< 

p 

O 

0 

u 

3 

0 

♦-  —  in 

*-  -  0 

—  in 

cr  0  0 

a  0  0 

cr  0  0 

UJ  0  0 

UJ  0  0 

UJ  0  0 

>  •  • 

>  *  - 

>  •  - 

d  0 

d  0 

d  0 

1  l 

d  d 

d  d 

0) 

l  1 

1  • 

Cn 

m 

a 

O  O 

0  0 

O  O 

p 

0  O 

O  0 

O  O 

O  O 

0  0 

O  O 

DO  — 

u  0  — 

DO  — 

VO  0 

v  0  0 

VO  0 

>  *  • 

>  •  • 

—  — 

>  •  • 

zoo 

zoo 

ZOO 

•y 

UJ  0  0 

Ui  O  O 

UJ  0  O 

♦-  0  0 

*-00 

♦-  0  0 

•  • 

< 

< 

< 

a 

z 

O  ui 

O  —  d 

»-  o 

a  2 

o 


Z  COMP  INCL  •  .000  THETA  •  .000  TALT  «  .0  RALT  ■  .0 

D1ST  amplituoe  phase  dist  amplitude  phase  dist  amplitude  PHASE  D 1  ST  amplitude  phase 

.000  104.7406  .0000  2.520  54.1079  62.0440  5.040  39.3953  64.96ie  7.560  33.4001  64.6874 


f-  —  — 

cd 

—  O 

in 

47 

o 

00 

V  CD 

to 

*-  o 

in 

CD 

Oi 

in 

© 

to 

to 

47 

TJ 

o 

cd  o-  oi  *-  in 

© 

in 

© 

CD 

Oi 

CD 

— 

o 

Oi 

•- 

CD 

© 

r- 

(O 

(O  Oi 

CD 

Oi 

00 

in 

Oi 

co  in  - 

o 

cd  o 

CD 

Or 

CD 

CD 

o  in 

Oi  O'  © 

to 

in 

•- 

o 

«- 

•- 

CD 

CD 

•» 

CD  in  to  *7  CD 

CD 

Oi 

o 

O 

© 

CO 

in 

CD 

o 

CD 

o* 

w 

CM 

<D  Oi 

CD 

Oi 

(D 

O 

47 

Oi 

© 

CD 

inn- 

©CD  ffl 

cd 

© 

T 

to 

in  cd 

o- 

O  CD 

o 

to 

to 

CD 

m 

47 

in 

o 

o> 

to 

soinoio 

O 

Oi 

m 

CD 

c 

O 

CD 

to 

in  to 

to 

CD 

*- 

<7 

r-  cm 

to 

Oi 

O' 

7 

O 

O- 

o  ©  O' 

o-  n 

o 

CO 

o- 

to 

to  to 

O' 

CD  O 

n 

in 

CD 

in 

<D 

m 

CD 

Oi 

o» 

Oi  oo  ©  ©  in 

— 

o* 

© 

CD 

to 

r>  cd  to 

©  © 

r- 

47 

Oi 

CD 

lO  47 

- 

CD  CD 

47 

Oi 

CT> 

O' 

in 

47  n  Oi  CM 

o  o  o 

o  o 

O  O  — 

Oi 

Oi 

Oi 

© 

© 

47 

c 

in  m  to  to  o» 

CD  CD  CD 

CD 

o 

Oi 

m 

V 

m 

to  to  O'  CD 

CD 

CD  O 

Oi  o*  r> 

* 

ID  ©  ©  © 

to  to  to 

CD  CD  tO  tO  CO  CD 

tO  CO  CO  CO  CO 

CO 

tO  CD  tO  CO  CO 

to  CD  (O  to  to  to 

tO  tO  CO 

to 

r» 

o- 

o* 

r- 

o- 

r- 

r- 

o* 

r> 

r- 

r-  t- 

r- 

o» 

CD 

CD 

CD 

CD 

CD 

0 

cd  cn  (\  o  tr  o 

^  in  MD  OUD  V 

47  <7 

©0'0©CD©©©©©  — 

Oi  CD  ©  47 

Oi 

<d  ©  © 

— 

Oi 

© 

Oi 

Oi 

©  Oi 

O 

_ 

m 

o 

© 

Oi 

Oi  ©  r> 

p) 

tDfOomor^ninnMD*-*' 

©  © 

—  O' 

r- 

oi  co  m  o  O'  ©  r- 

o  ©  © 

© 

© 

*- 

© 

Oi 

rv 

© 

o 

© 

Oi 

© 

O 

CO 

© 

© 

Oi 

© 

CJi 

C' 

CO 

©  n  © 

«— 

Oi  in  CD  O  Oi  <7 

ID  N  CD  CD  (D  ID  N 

©  © 

*- 

©  47 

o 

©*-©o©©cdo<©© 

© 

© 

O' 

© 

© 

© 

© 

Oi 

© 

O' 

<7 

»- 

m 

© 

<7 

© 

CD 

f" 

«- 

© 

©  — 

CO 

0 

to  ©  o  co  in  o* 

(j)  •-  n  in  n  oi  *• 

©  ©  o- 

©  O  Oi 

r>  ©  ©  ©  o>  o 

©  47  ©  © 

O' 

© 

© 

o 

- 

Oi 

© 

CO 

^7 

© 

© 

© 

O' 

©  © 

0- 

© 

o 

o 

o  — 

— 

— 

©  r>  rr  *7  47  47 

v  in  m  m  m  in  id  id  id 

© 

©  O' 

O' 

o-  o»  o-  o*  o-  ©  © 

©  ©  © 

©  © 

© 

© 

© 

© 

© 

© 

© 

©  © 

© 

© 

© 

© 

© 

©  © 

© 

o 

O 

o  o 

o 

o 

onnnoonnnnncjnnfonnn 

nnnnnnnnnnn 

© 

©  © 

© 

©  r> 

©  r> 

r> 

© 

©  © 

n  © 

ry 

CO 

CO 

n  co  *7 

<7 

* 

oooooooooooooooooooooooooooooooooooooooooooooooooooooo 

CDOniTMOfflOCilTMOCCOCMTnOfflOtN'IliJEOfMiJirCDOrtillDtDOri'ItOnOrM'JlDCCOrigiCCCOriiJ^OONV 

iTHTMOlDlOlDNt'r^NNODOCDtCCDaXJlCM^ClOOOOO*-*-  -  -  N  N  ^  (S  n  P)  O  f)  D  H  ^  47  <7  47  47  ©©©©©ID©© 


f«-  f'.  f««.  f*'  f'-  r"  0-  f'-  t-'-  0-  O'  0~  O'  0~  QD  CO  CD  CD  ©  ©  ©  ©  CD  ffl  0  CD  ©  CD  CD  ©  ©  CD  CD  ©©©©©©©©©©  CD  ©  ©  CO 


O' 

0©0©0i©©0i 

w- 

O 

© 

© 

© 

47 

O'  O' 

O 

© 

O 

© 

O 

© 

© 

O'  ©  47 

© 

© 

© 

© 

© 

© 

Oi 

47 

47 

O 

© 

O' 

© 

O' 

© 

© 

o 

© 

_ 

Oi 

n 

-  © 

*-conr)OicD'-o 

^7 

Oi 

O 

O' 

47 

© 

Oi 

Oi  © 

© 

Oi 

r- 

© 

© 

■*— 

O-  ©  © 

Oi 

© 

© 

© 

© 

© 

O 

Oi 

47 

© 

© 

© 

© 

•- 

© 

© 

© 

© 

© 

© 

47 

© 

Oi 

O'  © 

© 

fMNinTJTjlDOC 

O 

O' 

© 

© 

© 

© 

© 

O'  © 

© 

O' 

■*- 

© 

Oi 

O- 

47 

O  O-  47 

Ci 

© 

O' 

© 

© 

w- 

O 

© 

© 

47 

Oi 

o 

® 

© 

© 

O 

O- 

47 

o 

© 

Oi 

O' 

Oi 

© 

©  Oi 

<7 

o©4-o-r)©©oi 

© 

© 

Oi 

© 

© 

© 

O 

O'  47 

Oi 

© 

O' 

47 

Oi 

© 

O' 

©  Oi  O 

© 

© 

© 

© 

O' 

© 

Oi 

O 

© 

© 

47 

O 

© 

Oi 

O 

© 

© 

© 

o 

© 

© 

Oi  o 

© 

©©O'O'©©©© 

o 

_ 

Oi 

Oi 

© 

47 

©  ©  ©  O' 

O' 

© 

© 

O 

O 

•- 

Oi  ©  47 

47 

© 

© 

O' 

© 

© 

O 

_ 

Oi 

© 

47 

4J 

© 

© 

O' 

© 

© 

© 

O 

Oi 

r>  47 

©©©©©©©©O' 

O' 

o- 

O- 

o- 

O- 

O' 

f' 

O'  C" 

O' 

O' 

O' 

O' 

© 

© 

© 

©  ©  ©  © 

© 

© 

© 

© 

© 

m 

© 

CD 

© 

© 

© 

© 

© 

© 

© 

©  © 

© 

CD 

o 

o 

o 

O 

o  o 

© 

*-oinon©oo 

© 

O' 

47 

© 

47 

©  © 

O  O' 

© 

© 

CD 

O' 

© 

—  —  ©  © 

47 

© 

o 

47 

© 

© 

Oi 

O 

© 

© 

O' 

Oi 

47 

47 

•7 

© 

Oi 

© 

© 

Oi 

CD 

o 

Oi  © 

470i©©0—©0i 

© 

© 

© 

O 

O 

o- 

*-  o  © 

© 

© 

47 

o- 

© 

w~ 

4?  ©  © 

-- 

O 

© 

c 

X 

«7 

O' 

r- 

© 

o 

© 

© 

o 

o 

47 

© 

© 

© 

O' 

C 

© 

o 

o 

O' 

O  © 

Oi©CO©©OiO©0‘ 

47 

w- 

m 

© 

•— 

©  o#  o- 

w- 

©  © 

©  ©  ©  Oi  47  in  o 

©  O 

o 

o 

O 

O 

© 

© 

O' 

© 

© 

-- 

© 

© 

Oi 

© 

47 

o 

© 

o 

© 

o 

47 

0. 

—  r> 

©©0.©©0'-CiCi 

n 

47 

©  © 

©  C' 

O' 

©  © 

©  © 

© 

o 

o 

o 

— 

^  - 

Oi 

Oi 

Oi 

Oi 

Oi 

— 

— 

- 

— 

- 

O 

o 

o 

© 

©  © 

© 

© 

O' 

O' 

© 

© 

©  47 

© 

©©©©OOOOO 

©  O 

o 

o 

o 

o 

o  o 

o 

o 

o 

_ 

^  ^ 

_ 

_ 

_ 

_ 

_ 

o 

O 

o 

o 

o 

o 

o 

O 

o 

O  O 

nnconnTTi  ti  c 

47 

47 

47 

47 

47 

47 

47 

47  47 

47 

47 

47 

47 

47 

47 

47 

47  47  47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

47 

V  47 

O 

ooooooooo 

o  o 

o 

oooooooo 

O  O  O  O  O  O 

o 

o 

o 

o  o  o 

OOOOO 

o 

o 

o 

o 

o  o 

o  o 

o 

o  o 

o 

o 

o  o  o 

© 

© 

O  Oi 

47 

© 

© 

o 

Oi  47 

©  © 

o 

Oi  47 

© 

©  o 

Oi  4J 

© 

© 

O 

CN  47  ©  ffi 

o 

Oi 

47 

©  © 

o 

Oi 

47 

© 

© 

o 

Oi 

47 

© 

© 

o 

Oi 

47 

•o 

© 

o 

Oi 

<T 

© 

© 

o 

OI 

o  o 

4- 

4" 

— 

Oi 

Oi  Ci  Oi  Oi 

© 

©  © 

© 

©  47 

47  47 

47 

47 

© 

©  ©  ©  © 

©  © 

© 

©  © 

O' 

O' 

O' 

r- 

© 

© 

© 

© 

© 

© 

© 

G 

o. 

© 

o 

o 

o 

o 

o 

— 

4- 

© 

© 

5 

5 

© 

© 

© 

© 

5 

5 

©  © 

© 

©  © 

© 

5 

5 

5 

5 

© 

© 

© 

5 

5 

5 

5 

5. 

© 

© 

5 

5 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

©  © 

© 

© 

© 

w- 

Oi 

© 

© 

©  © 

© 

©  Oi 

Oi 

© 

Oi 

© 

47 

Oi 

Oi  — 

© 

Oi 

© 

© 

© 

47 

© 

© 

«- 

© 

o 

47 

47 

© 

Oi 

© 

o- 

47 

o 

© 

© 

47 

© 

r- 

Oi 

ffl 

© 

© 

© 

© 

© 

© 

© 

w- 

■»— 

CD 

© 

©  47 

47 

©  0" 

© 

•— 

© 

47 

© 

^7 

Oi 

©  Oi 

© 

Oi 

O' 

© 

O' 

•*- 

O 

47 

©  © 

Oi 

•- 

© 

© 

© 

© 

© 

© 

Oi 

47 

Oi 

© 

© 

© 

O 

C' 

47 

© 

© 

4J 

© 

Oi 

o 

© 

w~ 

© 

© 

Oi 

o  © 

© 

47  Oi 

o 

© 

0 

© 

© 

•— 

© 

©  47 

w- 

m 

47 

o 

© 

Oi 

V- 

© 

© 

© 

47 

47 

47 

Oi 

o 

O' 

Oi 

O' 

o 

Oi 

Oi 

— 

© 

47 

© 

© 

© 

© 

w 

47 

47 

V- 

© 

© 

© 

© 

© 

4- 

©  © 

47 

Oi  o 

© 

© 

© 

4“ 

© 

O* 

47 

Oi  O 

© 

© 

© 

© 

© 

© 

4" 

©  © 

© 

o 

O' 

47 

© 

© 

— 

© 

47 

-- 

O' 

© 

© 

47 

O 

© 

O 

© 

o 

© 

© 

© 

o 

© 

Oi 

n  47 

©  ©  ©  O' 

©  © 

©  o 

Oi 

Oi  ©  47 

©  ©  ©  t" 

© 

© 

©  o 

Oi  Oi 

©  47 

©  © 

©  O- 

O' 

© 

© 

© 

o 

_ 

Oi 

Oi 

© 

47 

47 

© 

5 

6 

© 

6 

7 

O' 

© 

© 

© 

©  ©  © 

©  ©  ©  ©  ©  © 

©  ©  © 

©  ©  © 

©  © 

©  ©  © 

©  ©  o  o  o  o 

O 

o 

o  o  o 

o 

o  o 

o 

o 

«- 

*- 

•- 

— 

-- 

— 

•- 

*- 

•- 

-- 

*- 

*- 

•- 

*— 

*- 

M— 

*- 

•- 

*— 

•— 

w- 

©O'©oi©47<D4y©oi©0'©  —  —  *-©f^©f^o4-‘©©©©©o©47^  ©©--©0'©oiO'Oi©©©©f^4T©4Tr>'^©©©4r 
0'^,©©©*'0i«-©©©4ro'©©©©0'C>i©0i0'©0'0i4r©©o©4T©©©©in©©©©ci©©o«-©o©  —  047©©© 
ooouiiN«D<rM(TMOnotooiiDn(i)r)or<(Dff)(NincooNnininiDiD(Dincro»-(T>NTj  --NnoiTjcDniDoroins  d  oi 
—  ooo©ooo>(D<DCD®©0'0'©©©©^©©pt.-»'-ocr>fflO'©©47©oi»-o©o-©©©oiOCDo-©47©  —  cdo*©© 


wwcwvnnnnr)r)nnnnr)nnnnroor*)f,)NryfNC><oir<fvf'<rvr<*-**T-»-^^^ooooooaicrioi(Ji 

©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©©^47^^ 


oooooooooooooooooooooo cvo oooooooooooooooooooooooooooooo 
VOlDOrt^lOCDOnVtOCDOOMJlOCDOnity^eDOWV^OOOW^tOCDOn^lOaDOCV^JDfflOfMTnOCDOWCtDCOO 
©©©iD©iD©©0'0'0-0'0'©©©©©cD(DrDCDcDooooo  —  —  *-^«p(r<nfnnnn«niiijij  TTCininminmc 


NwwwwwnwwNttMWNWPiwwPiMriMwnnnnnwnnnnnnnnnDnnnnonnnnnnnnop) 


ooi43o©©cT>©co)f'©©»-f,'0'OiCM©o©iD©0'©©0'*-©(n«'-o*-o©(N©0'470ior^(D©cD©<D^©©ooc*oi 
0'f,7r>o-©Oim©©*-©^,fnooiocr)(DO'0'iD4T©o©oi©©io*-0'"-o»’-©4jcD<DOiOiCr>  —  o©©©©0"®iD©to<Dc*> 
o^No>*-NN-o>Nfr)(D«intDUMroina)o>CD9cnooHfMO«n*“^c<inNnoDMD«-(Dcva)TjoincDO'(D^oO)0* 
v  —  a)ir>r>or>^ior^vor-r)aTm  —  r«-CNr»cNt^cNa)*-^T®'-<Tr^CD^rMr>nrotNoa)inrMNcNincDooC7>r-^Tovr'Oi 

c*r>r>vin©t0O'®®(D©o*-  —  c*nr>v47©©©<0o*t^i^©®®®tD<ncD(D<D(DCD®©©O'O'<D©©vn—ocDo-ino 


f^o*-n©~  —  ©•»-in<0(D4ytoor*ix>iDiDcnr>0'47t-®o®0'Oiotr>0'OiO'Oi©to©oi*-<T>CNnoi©©tD^©c*©  —  ®o 
tn<D®©©©oir>o©oio©  —  ©©0'Orio#ffl<D©4rto~©tn  —  <D©no©Ofn©iD©©r>  —  ©©©^©©©©ci  —  vm 
niocnoino'n^diow^i^in^  ©  o~  •-  ©  o  ©  —  ©i/iwoMnnoaiiOii  —  cd  ©  m  o  47  *-  ©47  »-a)iono(oo»-n(D»- 
r«*tO(8ifl^<o©«tf>oift  -  <D<M(DTror**r>04Dr)o>to<*>oU50or^r>Of^R3’or^^--r^v*-r^^T*-f*‘<s,^CDinri((T>tow*“ 

r^ffnDM£)ininwnr)op(^^^ooo(j'tJ'coa5COffiNNN^iciOifMninirc«Tnr)nrtrtN»'*-'-oooa>(T<o>0»  i 

^NaMOtfWtfXflOMD^UMDiDwajiocoiiKcinininininin^jrinininninin^ininintnjft^intfvini/vioi/Jinintnvvvv 


oooooooooooooooooooooooooooooooooooooooooooooooooooooo 
nviocDOciviocDonvtfeDort^tDQorwtoisoncicoon^iDoo^wcDortnioaion^uMDONvoQ 
oooo»-*>*-«-*-rtMP*wr<i,l)r)nnr)ij^  47  v  47©©©©©iDH>ir>©iD0'f^r‘'0'  o»  mcDixjcDCDOitT'cnffia'OOooo 


•4  ^ 


Figure  5.  Second  page  of  output  produced  by  lines  1-24  of  sample  deck 


1 


segment.  The  headings  "VERTICAL",  "BROADSIDE",  and  "END  FIRE"  are  for  the 
magnitude  in  dB  and  phase  in  radians  of  the  components  of  the  excitation 
factors.  The  headings  "V",  "B",  and  "E"  are  for  the  relative  magnitude  of 

these  components  at  the  beginning  of  each  slab.  This  list  appears  only  when 
ICOMP,  TALT,  RALT,  or  the  mode  parameters  are  changed.  In  figure  5  we  see  the 
list  of  computed  amplitude  in  dB  above  1  )V/m  and  phase  relative  to  the  speed 
of  light  in  degrees  printed  as  functions  of  distance  in  Mm.  This  list  is 
preceded  by  a  line  of  print  indicating  the  received  field  component,  the 
transmitting  antenna  orientation,  and  the  altitudes  of  the  transmitter  and 
receiver.  The  plots  for  these  calculations  are  shown  in  figure  6. 

Lines  25-31  of  figure  3  call  for  the  Ey  fields  and  phase  plots.  Three 
transmitting  antenna  orientations  are  specified.  The  resulting  plots  of  the 
amplitude  and  relative  phase  are  shown  in  figures  7a  and  7b,  respectively. 

There  are  apparently  only  two  curves  in  figure  7a  because  the  amplitude  of  Ey 
from  a  vertical  source  is  the  same  as  that  from  a  horizontal  end-fire  source. 

Lines  32-39  of  figure  3  call  for  a  rotating  transmitter  antenna  inclined 
at  45°.  The  calculations  are  to  be  made  for  two  distances.  The  first  page  of 
printout  generated  by  these  lines  is  shown  in  figure  8.  The  mean  signal  and 
the  standard  deviation  of  the  signal  for  the  rotation  are  printed  at  each  dis¬ 
tance.  The  table  of  amplitude  and  relative  phase  as  a  function  of  rotation 
angle,  R,  appears  as  a  result  of  NPRINT  =  1.  The  plots  resulting  from  the 
calculations  are  shown  in  figure  9. 


Lines  40-74  of  figure  3  illustrate  the  IOPT  =  4  option, 
resulting  from  these  calculations  are  shown  in  figure  10. 


The  plots 
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Figure  7b.  Plotted  relative  phase  output  generated  by  lines  25-31  of  the 
sample  deck 
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Figure  9.  Plotted  output 
generated  by  lines  32-39  of 
the  sample  deck 
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CALL  SYM80L(XP, VP, . 1 ,0T,0. , 16) 

XP«0. 

YP=VP-. 15 

1MNPL0T-2)  101,102,103 
101  CALL  SYMBOL! XP.VP, .1 , 9HAMPLI TUOE , 0 . ,9) 


UL ( 1 )«. FALSE. 

UL(2)». FALSE. 

CALL  CURVE  (XL,  YL.Ijl,  2,0. .0. ,1. ,1 . ,KK) 
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APPLIED  PHYSICS  LABORATORY 
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